Quantum criticality in inter-band superconductors 



O 

(N 

m 

(N 



Aline Ramires, Mucio A. Continentino 

Centro Brasileiro de Pesquisas Fisicas, 
Rua Dr. Xavier Sigaud, 150 - Urea, 
22290-180, Rio de Janeiro, RJ, BraziQ 

In fermionic systems with different types of quasi-particles, attractive interactions can give rise to 
exotic superconducting states, as pair density wave (PDW) superconductivity and breached pairing. 
In the last years the search for these new types of ground states in cold atom and in metallic systems 
has been intense. In the case of metals the different quasi-particles may be the up and down spin 
bands in an external magnetic field or bands arising from distinct atomic orbitals that coexist at 
a common Fermi surface. These systems present a complex phase diagram as a function of the 
difference between the Fermi wave- vectors of the different bands. This can be controlled by external 
means, varying the density in the two-component cold atom system or, in a metal, by applying an 
external magnetic field or pressure. Here we study the zero temperature instability of the normal 
system as the Fermi wave-vectors mismatch of the quasi-particles (bands) is reduced and find a 
second order quantum phase transition to a PDW superconducting state. From the nature of the 
quantum critical fluctuations close to the superconducting quantum critical point (SQCP), we obtain 
its dynamic critical exponent. It turns out to be z = 2 and this allows to fully characterize the 
SQCP for dimensions d > 2. 
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In strongly correlated materials superconductivity can 
be suppressed in different ways. Most commonly, this 
is accomplished by an external magnetic field, applied 
pressure or dopingir— . However, the point in the phase 
diagram where the critical temperature T c vanishes as 
a function of the external parameters is not necessarily 
associated with a SQCP. In the case of superconductiv- 
ity induced by antiferromagnetic fluctuations due to the 
proximity of an antiferromagnetic quantum critical point 
(AFQCP) 8 , as the system moves away from the AFQCP, 
these fluctuations change from attractive to repulsive and 
superconductivity just fades away^. 

Inter-band superconductivity is a long standing prob- 
lem in condensed matter physic o 10 ' 11 . It can be real- 
ized, for example, by applying an external magnetic field 
in a metallic system. The field splits the band and the 
phase diagram as a function of the mismatch of the Fermi 
wave- vectors of the up and down spin bands can be in- 
vestigated. It also occurs for the case of superconduc- 
tivity mediated by valence fluctuations, where the most 
relevant correlations are due to inter-band pairing^. 
In multi-band metals, pressure can be used to change 
the hybridization and the mismatch of the Fermi wave- 
vectors of the different bands coexisting at the common 
Fermi surface 1 ^. More recently this problem has received 
much attention due to the possibility of investigating it 
in cold atom systems. In this case the attractive interac- 
tion between two fermionic systems with different Fermi 
wave- vectors can be tuned by Feshbach resonance^ 4 -. The 
Fermi wave- vector mismatch Skp in this case is con- 
trolled by varying the density of the atoms and this al- 
lows the phase diagram to be fully explored. In the core 
of neutron stars, up an down quarks, in different num- 
bers, with attractive interactions may give rise to color 
inter-band superconductivit y 13 ' 14 . We neglect here or- 
bital effects. Whenever a magnetic field is mentioned, it 
is to be considered as an external agent whose only effect 



is to produce the mismatch of the Fermi wave-vectors. 
While avoiding complications related to the vortices 1 ^, 
the present approach is relevant for many interesting 
problems^ 4 - as pointed out above. 

The zero temperature phase diagram for inter-band su- 
perconductors with s-wave pairing has been established 
using mean-field calculations 11 . As the Fermi wave- 
vector mismatch increases, there is a first order phase 
transition from a BCS-like state to an FFLO or PDW 
state. As the mismatch further increases there is a con- 
tinuous transition to a normal metallic phase* 1 -. As new 
systems are discovered and the possibility of finding ex- 
perimentally new exotic PDW phases increases, it is im- 
portant to fully understand the nature of the normal to 
PDW transition beyond the mean-field approximation. 
Here we investigate this transition using a new approach. 
It is based on a perturbation theory for retarded and ad- 
vanced Green's functions 1 ^. We relate the appearance of 
superconductivity in the multi-band system to the diver- 
gence of a generalized susceptibility^ x{li w ): like in the 
Stoner criterion for ferromagnetism. 

We start in the normal phase and calculate the re- 
sponse of the system to a frequency and wave-vector 
dependent fictitious external field which couples to the 
superconducting order parameter—. For simplicity we 
consider here the case of s-wave superconductivity. 

In our approach, at the level of the RPA approximation 
presented here, we have to calculate only single-particle 
Green's functions. This is different and simpler than the 
usual linear response theory, which relates the response 
of the system to two-particles Green's functions 1 ^. 

At zero temperature, as the Fermi wave-vectors mis- 
match is reduced, we find a divergence of the static part 
of the susceptibility at a finite wave- vector q — qo- This 
occurs at the critical mean-field value of the mismatch 
which destroys pair density wave superconductivity. This 
divergence implies that even in zero field the system can 



2 



have a finite inhomogeneous superconducting order pa- 
rameter. The condition for the superconducting instabil- 
ity can be expressed in the form of a Stoner-like criterion, 
1 — Uxo(q = Qo,Skp) = 0. This determines either a crit- 
ical value of the attractive interaction U c above which 
the system is superconductor or a critical mismatch be- 
low which superconductivity sets in. The coupling of 
the order parameter to the electronic degrees of freedom 
determines the frequency dependence of the generalized 
susceptibility and the nature of the quantum fluctua- 
tions, i.e., the dynamic critical exponent at the PDW 
SQCP. The theory is equivalent to a quantum Gaussian 
approach and since the dynamic exponent turns out to 
be z = 2, it yields the correct description of the SQCP 
for dimensions d > 2. 

We start with the following Hamiltonian describing a 
two-band system with inter-band attractive interactions, 



(i) 



For simplicity, we omitted spin indexes 14 . The bands a 
and b can be the up and down spin bands of a metal polar- 
ized by an external magnetic field ftii, different types of 
atoms or the hybridized bands of a multi-band metal^. 
The interaction U is an attractive interaction between 
the fermions in different bands. We calculate the re- 
sponse of the normal two-band system to a wave-vector 
and frequency dependent fictitious field that couples to 
the superconducting order parameter of interest, in this 
case, 
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where the frequency ujq has a small positive imaginary 
part to guarantee the adiabatic switching on of the field. 
The response of the system to the fictitious field, will 
be obtained using perturbation theory for the retarded 
and advanced Green's functions^. We start in the nor- 
mal phase where the superconducting order parameter is 
zero in the absence of the external field g. We split the 
Green's functions, normal and anomalous, in two contri- 
butions. The first of order zero and the second of first 
order in the field g. For the anomalous Green's functions 
<C aLl^jo-' f° r example, we write, 
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where e£ = ^jt^e^-^i ri > , (a = a,b). For the first 
order Green's functions we get, 
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The first three first order propagators are written in 
terms of zero order propagators associated with super- 
conductivity. Since we are analyzing the system in the 
normal state, these anomalous zero order propagators 
vanish and consequently the first order ones also van- 
ish. The three last first order propagators are written 
in terms of conventional zero order propagators, but just 
the last one is different from zero. It is given by, 



«4i4'»L 1): 



>k+q,k' 



{U5A? + g) 



2tt( W + e£)(w + uj - e b k+q ) 



(4) 

We use the fluctuation-dissipation theorem^ to obtain 
the correlation function <5A°° from the above anomalous 
Green's functions, we get, 
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In the normal phase and in the absence of the fictitious 
field, the relevant zero order Green's functions can be 
easily calculated, 
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where F u {G(u>)} = -J dujf(uj)[G{uj + ie) - G(uj - ie)] is 
the statistical average of the discontinuity of the Green's 
functions G(w) on the real axis 1 ^. The function f(oj) is 
the Fermi-Dirac distribution. 

We define the generalized inter-band susceptibility by, 
X a o b (q,uj) = 2irj: k F u {Gt a (-u J )G k b +q (o J + u Jo )} or, 
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The superconducting response to the external field is ob- 
tained in the form, 



SA ab 



(6) 



Finally, the susceptibility Xq & ((?,u;), (Eq. O, is given by, 



2tt — el_ q + t» k+q - uj q 



(7) 



Let us consider the case of a metal polarized by an 
external magnetic field h. In this case the a and 6-bands 
are given by the dispersion relation of the up and down 
spin bands, respectively, 
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where kp is the original Fermi wave-vector of the unpo- 
larized band (h — 0). The mismatch of the new Fermi 
wave- vectors is given by, 5k f — k F — k h F — h/vF, where 
vf is the Fermi velocity. At zero temperature, the gen- 
eralized susceptibility Xo b (2?, = 0) can be calculated 
and we get, 



X 8 b (2 g ,0)=p{l-ln(^A_)-I 
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where q = VFq/h = q/(k F - k b F ) and p = 3/(8tt 3 Ef) 
is the density of states of the unpolarized system. The 
quantum critical point separating the normal metal from 
the superconducting PDW phase is obtained from the 
condition, 1 — Uxo b (2q,0) = 0, which can be written as, 
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where Ao = vpkc exp(— 1/pU). Notice that a PDW state 
is only possible for q > 1. The equation above determines 
the critical field (or mismatch) for a fixed value of q. 
This is given by, 
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The maximum value of /ijv for which the instability 
occurs is denoted by he- It is easily obtained by differ- 
entiating h jv (q) with respect to q. From the equation 
(dlnliN / dq)q = q c — 0, we find that: 
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q c = -ln(^L 
° 2 \q c 

This gives q c = 1.2, that substituted in Eq. [10] for hN 
yields he = h^iqc) — 1-5 Aq. This agrees with the result 
of Takada and Izuyama— for the vanishing of the FFLO 
phase starting from this phase. 
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FIG. 1: Schematic T c x h phase diagram for the inter-band 
superconductor. At zero temperature, as the external mag- 
netic field or Fermi wave-vectors mismatch is reduced (indi- 
cated by the arrow), the normal metal becomes unstable to 
an inhomogeneous superconducting phase. This is signaled by 
the divergence of the g-dependent generalized susceptibility. 
The FFLO superconductor quantum critical point at he has 
a dynamic exponent z — 2 and this allows to determine the 
behavior of the thermodynamic quantities along the critical 
trajectory (CT) as given in the text. As the external field 
is further reduced there is a first order transition from the 
FFLO state to a BCS superconductor^ which is not studied 
here. The crossover line T x marks the onset of Fermi liquid 
behavior with decreasing temperature. 



As discussed in RefJ^ an FFLO phase can also be 
induced by hybridization (V). The results above can 
be easily extended for the case of hybrid bands 13 . The 
Fermi wave-vectors mismatch is given by, k F — k F = 
4V/vfO- + a )j where a < 1 is the ratio of the effective 
masses of the quasi-particles. The mismatch is propor- 
tional to the hybridization V (or pressure), which now 
plays the role of the magnetic field. 

In order to study the fluctuations close to the PDW 
SQCP we have to expand Eq. [7] for small frequencies ujq, 
h = he and q = qc- At the same level of approxima- 
tion for q of Eq. [SI we get for the denominator of the 
generalized susceptibility in Eq. [6[ 
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The coupling of the superconductor order parameter to 
the electronic degrees of freedom gives rise to Landau 
damping of the superconductor quantum fluctuations. 
These modes are purely evanescent with an imaginary 
dispersion relation. 

As in Hertz approach 20 to quantum phase transitions, 
we can construct from the dynamical susceptibility a 
quantum Gaussian action for this problem 2 ^, which can 
be written as, 



S= d d q duj[g + uj + q 2 ] \ip(q,u)\ 



(12) 



where g — h — he measures the distance to the PDW 
SQCP, tp(q, lj) is the PDW superconductor order param- 



4 



eter and q — qc — > q. The scaling properties of the free 
energy associated with this action allows to identify the 
dynamic exponent of the SQCP at which the FFLO in- 
stability occurs. This is easily found as z — 2. The 
knowledge of the dynamic exponent and the quantum hy- 
perscaling relatio n 21 ' 23 , 2 — a = v{d + z), allows to fully 
characterize the universality class of the FFLO quantum 
phase transition. Since d e ff — d + z, for d = 2 and 
d = 3, the critical exponents are Gaussian with possible 
logarithmic corrections for d = 2. The reason is that for 
d > 2, interaction between the fluctuations are irrelevant 
in the renormalization group sense and the Gaussian ac- 
tion, Eq. 1121 gives the correct description of the quantum 
phase transition^. 

Next step in the calculations is to extend the results 
to finite temperatures and particularly to obtain the cor- 
rection to the mean-field value of the shift exponent ip 
of the critical line at very low temperatures close to the 
PDW SQCP. This can be carried out in different ways^ 4 -, 
such that, T c (h) oc \h - h c \ 1/,p . It turns out that 24 
ip = z/(d+z—2) for d e ff > 4 and can be determined from 
our knowledge of the dynamic exponent, which yields 
ip = 2/d. In the same way, scaling^ gives the contribu- 
tion of quantum critical fluctuations to the specific heat, 
C/T oc T^-^/ 2 , along the critical trajectory (CT) shown 
in Fig. [TJ Also shown in this figure is the crossover line^i, 
T x oc \h — hc\, below which the system enters a Fermi 
liquid regime. We have used that vz = 1 for d > 2. 

We have studied the zero temperature instability of a 
normal two-band system with inter-band attractive inter- 



actions as the Fermi wave-vectors mismatch is reduced. 
The model describes a superconductor in an external 
field in the absence of orbital effects, a cold atom system 
with two atomic species or a two-band metal with mixing 
tuned by pressure. We find that as the mismatch of the 
Fermi wave-vectors is reduced, the system becomes un- 
stable to an inhomogeneous, FFLO-like, superconducting 
phase characterized by a wave- vector q = q c . We have 
introduced a new method to deal with systems coupled to 
a space and time dependent external field and respond 
adiabatically to this perturbation. The appearance of 
superconductivity is related to the divergence of a static 
generalized susceptibility for q = q c at a SQCP. The 
results obtained are valid to first order in the perturba- 
tion and coincide with those of linear response. How- 
ever, in our approach the basic elements to be calculated 
are single particle Green's functions and not the usual 
two particle propagators of linear response theory. The 
present theory extends the mean-field treatment to in- 
clude fluctuations close to this SQCP. We have obtained 
the dynamical critical exponent at the PDW SQCP and 
from that we have identified the universality class of this 
quantum critical point for d > 2. This allows several pre- 
dictions to be made for the thermodynamic and transport 
properties close to the quantum phase transition and for 
the shape of the critical line for dimensions d > 2. 
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